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ABSTRACT A generalisation of the classical Calogero-Moser model obtained 

(^ • by coupling it to the Gaudin model is considered. The recently found classical 

"^ ! dynamical r-matrix [E. Billey, J. Avan and O. Babelon, PAR LPTHE 93-55] for 

the Euler-Calogero-Moser model is used to separate variables for this generalised 

Calogero-Moser model in the case in which there are two Calogero-Moser par- 

A ' tides. The model is then canonically quantised and the same classical r-matrix 

D ! is employed to separate variables in the Schrodinger equations. 

> 

X 

■^ ■ 1 Introduction 

Separation of variables in the Hamilton- Jacobi equation is one of the methods of solving 
completely integrable models of classical mechanics (see e.g. [Q). If the model admits a 
Lax pair representation and hence is described by a classical Yang-Baxter algebra 0, the 
separation of variables can be achieved by the functional Bethe Ansatz 1|^ . By a classical 
Yang-Baxter algebra we mean an infinite-dimensional Poisson bracket algebra given by an 
N X N matrix L{X), of which components are dynamical variables, and an A^^ x A^^ matrix 
ri2(A,/i), which may depend on dynamical variables [§] [0 01 and is such that 

{L^iX),L2i^^)} = [n^lA,/.), Li(A)] - [r2i(/i, A), L^l/i)]. (1) 
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Here X, fi are free parameters, -^i(A) = -^(A) / and L2{fi) = / ® L{fi), and r2i(A,/i) = 
Pri2(A,/i)P, where P is an A^^ x A^^ permutation matrix. L{X) is called a Lax matrix 
and ri2{X,fi) is called an r-matrix. Matrix ri2(A,/i) has to satisfy a relation necessary to 
ensure that the Poisson bracket (|I]) obey the Jacobi identity. This relation can take a form 
of either the classical Yang-Baxter equation (see e.g. []TU|) or its dynamical generalisations 



lUlllEIl if '^i2(A,/i) depends on dynamical variables. We deal with this last case in the 
present paper. 

Once the classical dynamical system is written in the form (0) , we can define a generating 
function for integrals of motion, 

t(A) = hrL'iX). (2) 

From (|I]) it follows that {t(A),t(/i)} = 0, hence the integrals generated by t(A) are in 
involution and the system is integrable. The main goal of the functional Bethe Ansatz is 
to separate the variables in the Hamilton- Jacobi equations for the integrals generated by 
t(A). 

For models with sl{2) symmetry the functional Bethe Ansatz can be described as follows. 
We begin with a Lax matrix 

/A(A) C(A) \ 
^^^^-U(A) -A{X)) ^^' 



and, for simplicity, we assume the following form of a classical r-matrix 
ri2iX,n) 



fa{X,fi) \ 

6(A,/i) 

c{X,fi) 

V a{X,fx)J 



(4) 



for some functions a(A,/i), 6(A,/i) and c(A,/i). Functions a(A,/i), 6(A,/i) and c{X, fi) satisfy 
relations necessary to ensure that the Poisson bracket (|l]) obey the Jacobi identity. The r- 
matrices of this type appear in several integrable models such as homogeneous spin chains, 
heavy tops and systems on Riemannian manifolds of constant curvature. For any Lax 
matrix (|^) and any r-matrix (^ the Poisson bracket algebra (^ takes the form 

{A(A),S(/i)} = 6(A,/x)S(A) + a(/i,A)S(/i), 

{AiX),Cifi)} = -c(A,/i)C(A)-a(/i,A)C(/i), (5) 

{5(A), C(/i)} = 2c(A,/i)A(A) + 26(/i,A)A(/x), 



and zero for the remaining brackets. The separation variables are defined by 

B{X,) = 0, P, = A{Xi). (6) 

Using the algebra (||) we find [Q 



{X„ X,} = {P„ P,} = 0, {X„ P,} = 5,, ^lim^ ""^^^xT^ ^ 

so that a sufficient condition for (Xj, Pj) to be canonical variables reads 

a(A,X,)P(A) 

hm ; : = 1, 

A-x. B'{X,) 
or equivalently 

a(X, + /i,X,) = -^ + 0(/i). (7) 

We stress that condition (|^) remains the same for systems described by non-dynamical 
and dynamical r-matrices. This in particular allows one to apply the functional Bethe 
Ansatz procedure to the latter, as discovered in [|Il|]. This observation seems to be of great 



importance and interest, since there are integrable models, such as the Calogero-Moser 
model, for which the non-dynamical r-matrices do not exist. The models of this kind have 
recently attracted much attention and the theory of dynamical r-matrices requires closer 
investigation [Q. 

Having defined the separation variables we easily find the separated equations 

P^ = t{X,), (8) 

which follow directly from (^. 

In this paper we employ this procedure of separation of variables to solve the generalised 
Calogero-Moser model. In order to do so we use the recently introduced r-matrix |^ which 
depends on dynamical variables. In the next section we describe a generalisation of the 
Calogero-Moser model which is achieved by coupling the Gaudin model to it. Then we 
specialise to the case in which there are two Calogero-Moser particles. In this case the 
system has sl{2) symmetry and is governed by an r-matrix of the form @ which in addition 
satisfies condition (|^), so that the separation of variables procedure can be apphed. We 
do it in Section 4. In Section 5 we quantise the model and use the quantum version of 
the functional Bethe Ansatz to separate variables in the Schrodinger equation [0. The 
quantum functional Bethe Ansatz we use is obtained by the canonical quantisation of the 
classical procedure just described. 



2 Description of the Model 

In Billey, Avan and Babelon proposed the parameter dependent classical dynamical 
r-matrix for the generalisation of the Calogero-Moser model constructed by Gibbons and 
Hermsen [0]. The model is governed by the Hamiltonian 

where the dynamical variables (gi,Pi)i=i,...,Af, {fij)i,j=i,...,N satisfy the Poisson bracket alge- 
bra 

Ui, Qj} = {Pi,Pj} = 0, {qi,Pj} = 5ij, {fij, fki} = Sjkfii - Sufkj- (10) 

The model is integrable, when restricted to the surfaces {fa = const)j=i^...^Ar. In this section 
we propose a further generalisation of (P) which is achieved by introducing a non-trivial 
internal structure to the dynamical variables {fij)i,j=i,...,N- We consider a system governed 
by the Hamiltonian (Rf) but with the dynamical variables {qi,pi)i=i n, (/iw)»i=i n, which 

' ' ■' a=l,...,M 

satisfy 

{ft,JM} = Sap{S,,fti-Saf^,) (11) 

and all but the last of relations ([T0|) . The dynamical variables {fij)ij=i,...,N are now defined 
by fij = J2a=iftj- To see that the model (|)-(0) is integrable when {fa = const) 
we construct the corresponding Lax matrix 



i=l,...,N 



N / M fa \ N / f M fa \ 

^(A)=EU + Et^ e.+ E -^ + Et^ e,, (12) 

i=l\ a=l'^-^aj r,j=i \(li - (Ij a=l-^'-^o.J 

where {eij)ki = ^ik^ji and ei > €2 > ■ ■ ■ > ^m are arbitrary parameters. This form of 
L(A) immediately reveals that the generalisation of the Calogero-Moser model we discuss 



here is achieved by the coupling of the M-particle Gaudin model [15| to the Calogero-Moser 



system. Using equations (|T0|) and (ITT]) we can write the Poisson brackets of the components 



of the Lax matrix (|T2|) in the unified form 



N f.. _ f.. 

{Li(A),L2(/i)} = [ri2(A,^),Li(A) + L2(/i)] - E /"_ " e,,^eji. (13) 

The r-matrix here is the dynamical one considered in [§] 

ri2(A, Ai) = E 1^ ^ ^«i ® ^i» + 1^ E ^«« ® ^««- (^^) 



We note that ri2(A,/i) = — r2i(/i, A). From ( ]T3| ) we immediately learn that the model (||) 
is integrable provided fa = fjj, i,j = 1,...,N, since in this case the Poisson bracket 
algebra (|I3]) takes the form (|1|). In particular, the model is integrable on the surfaces 
{fa = const)j=i AT. This final reduction is possible, because each of fa Poisson commutes 
with the Hamiltonian H. The Hamiltonian H is recovered as H = ^ Jq 2^trL^(A), where 
A is considered as a complex variable and the contour C encloses the origin. 

3 The N=2 Case 

Now we focus on the N = 2 case for which we write a complete set of integrals of motion. In 
the next section we will use the functional Bethe Ansatz procedure to separate the variables 
in the Hamilton-Jacobi equations corresponding to these integrals of motion. 

In the centre of mass frame, the Lax matrix in this case has the form (^), with 

A/ qa C M qa q M qa 

A(A)=p+Et^, 5(a) = -^ + E/-^ ^(a) = ^ + E/V ^^^) 

a=l ^ ~ ^a Q a=l ^ ~ ^a 9 a=\ ^ ~ ^a 

where 5^, 5^ are generators of the so(2, 1) Poisson algebra, 

{S3", si} = ±5^pSL {SI, 5^} = 25^pSi 
and {q,p) are the relative coordinates, {q,p} = I, q = qi — q2, p = \{pi — ^2)- The matrix 



([T5| ) is obtained from (|T2|) by subtracting the centre of mass motion. The non-zero Poisson 
brackets of the components of the Lax matrix (|T5|) take the explicit form 

5(A) -5(/i) B{X) 



{A{X),B{^^)} 



X — fi q 



{AA,C/i} = + , 16 

X — jj, q 



{5(A), c(/i)} = 2 ( ^'""i ' ^'^> - ^^"^ : ^'^> + ^ 



'A{ X)-A{^i) _ A{X) - A(/i) ^ S,' 
A — /i q q^ 

The Poisson algebra (16) can be written in compact form (|l|) with ri2(A,/i) of the form 



where 

«(-^> /^) = T > KK Ai) = T \- -, c(A, /i) = , 

A — /i X — ft q X — fi q 



provided that 6*3 = 0. This follows immediately from ([T3|). Writing explicitly the generating 
function (|^), 

,,-, 2_S-S+ ^ ^3-^3^ + SZSJ , ^ 2pS^ + l(g!!g+ - S'^S.) 



5^ a^fcl (^ - e«) (^ - e/3) ^1 -^ - ea 



we can compute that {t(A),S'3} = 0, hence the reduction S3 = can be performed. A 
complete set of integrals of motion H, Ga, H^, a = 1, . . . , M is found by rewriting the 
generating function t(A) in the following form 

M TT M (^ 



a=l 



^ ~ ^a a=l (-^ ~ ^a) 



where 



2 S-S-i 
H = p 



G„ = [S^f + S^Sl, (18) 

M o qa ql3 , ( na qf^ < qP qa \ i 

13=1 e„ - 6/3 q 

Integrals of motion (^) are in involution provided that S3 = 0. Notice that if this is the 
case, then J2a=i ^a = 0. The first of the integrals of motion ( [Tsl) is the Hamiltonian of our 
system and for each a = 1, . . . , M, Ga is a Casimir function hence it remains constant on 
each symplectic leaf of the manifold on which the system is realised. Therefore we have 
M independent integrals of motion. We can represent the variables S" in terms of the 
canonical coordinates and momenta {xa,Pa), {xayPp} = Sai3, a, /3 = 1, . . . , M as follows 

Ss=l^aPa, SI = ^pI S^^ = ~xl (19) 

In the representation (^) the first integrals ([I8D take the form 

d2 m 

Ga = 0, a = l,...,M, (20) 

1 ^^ M^ 1 ^ 

Ha = -7 E — ^^ + p^c^Pa + :r E (pl^i - ^IpD^ 

where R^ = Yla=i ^a ^^nd Map = PaXp — XaPp- The constraint 6*3 = translates to 

M 
E ^c.Pa = 0. (21) 

If we choose H as the Hamiltonian of the system, as in fact we are doing here, then the 
constraint (|21|) implies that R^ = const and we have a two particle Calogero-Moser model 
coupled to the free motion on the sphere S'^~^. 

6 



4 Separation of Variables 

The dynamical system (120) can be solved by separation of variables in the Hamilton- Jacobi 
equations. The separation is carried out in the framework of the functional Bethe Ansatz 
IP] as described in the introduction. First we consider the first of Eqs.(|^), which in our case 
reads 

q i^i A - e„ 
It is a polynomial equation of degree M and it has M different solutions Xi, i = 1, . . . , M. 
We use the Vieta theorem to derive the following expressions for q and x^, a = 1, . . . , M 

A/ M p2 ttM ( y _ ^ \ 

i=i Q=i q 11/3=1 ie/3 - eaj 

We can choose the roots Xj, -i = 1, . . . , M in such a way that Xi > . . . > Xm- From ( P^ ) 
we then learn that the separation coordinates Xi, i = 1, . . . ,M satisfy the inequalities 

Xm < ca/ < Xm^i < . . . < Xi < ei or cm < ^m < . . . < ei < Xi. (23) 

Since a(A, /i) obeys (|^ it is clear that the canonical momenta Pi, i = 1, . . . ,M can be 
defined by the second of Eqs.(^. Therefore Eqs.(^) give a full set of canonical variables 
(Xj, Pj)j=i,...,M 5 {Xi, Pj} = 5ij, and we can proceed to separation of variables. We seek the 
common solution S{Xi, . . . , Xm) of the Hamilton- Jacobi equations 

r(^^ ^^ X x\-F H (^^ ^^ X x\-F 

in the separated form S{Xi, . . . , Xm) = Y^fLi Si{Xi). Here Ea are such that Y.a=i ^a = 0. 
From (H) we find the separated equations 

(t)^-^-l^-O, .^1,...,M. ,24) 

If we consider only a flow generated by the Hamiltonian H, then E^, a = 1, . . . ,M have 
the meaning of arbitrary separation constants. 

The separated Hamilton- Jacobi equations (0) can be also used to express the Hamil- 
tonian H in terms of the separation coordinates (Xj, Pj)j=i ... jvi- Eliminating constants E^, 
a = 1, . . . ,M from Eqs.(p^), we find that 

Ei=l Pi I\a=liXi - ea) l\j=i{Xi - Xj)~ 



Ei=l Ila=liXi - ea) Ilj=iiXi - Xj] 



Finally we would like to stress that to derive separated equations such as Eqs. (p^ ) we 
do not have to specify to the representation (|1^) of S", a = 1, 
can also separate the Hamilton- Jacobi equations for integrals (^ 
separated equations read 

2 M 



dSi 



E^ 



Gr 



dX, 



0=1 



X,; 



{Xi 



, M. In other words we 
In this general case the 

:0. 



5 Quantisation 



The classical system (j^) can be quantised in the canonical way by replacing the Poisson 
brackets { , } with the commutators —i[ , ]. Also, the functional Bethe Ansatz method 
can be quantised in this way [1^ . We can consider the Lax matrix (^ which generates the 



Gaudin algebra (see e.g. [p!6| ) 



[Li(A), L2{fi)] = i[ri2(A, /i), Li{\)] - i[r2i(/i, A), L2(/i)]. 



(25) 



A generating function t(A) for commuting integrals of motion is given by (0). We consider 
r-matrices of the form (^). Functions a(A,/i), 6(A,/i), c(A,/i) may depend on dynamical 
variables, hence they are operator valued functions in general. They satisfy similar condi- 
tions as before to ensure that the commutator (^) satisfies the Jacobi identity. In contrast 
to the classical case however, functions a{X,fi), 6(A,/i) and c(A,yu) must satisfy additional 
conditions to make the Gaudin algebra (p5|) consistent. The sufficient ones are 

[a(A,/i),A(A)] = [a(A,/i),i?(A)] = [a(A,/i), C(A)] = 0, 

[6(A,^),5(A)] = [c(A,/i),C(A)]=0. (26) 

If the conditions (p6|) are satisfied, then the algebra (|25| ) takes the explicit form 

[A(A),5(^)] = ?(6(A,^)5(A) + a(/i,A)i?(/i)), 

[A(A),C(/i)] = -z(c(A,^)C(A) + a(/i,A)C(^)), (27) 

[5(A), C(/i)] = ^([c(A,^),A(A)]+ + [%,A),A(/.)g, 

where [ , ]+ denotes the anticommutator. 

The separation of variables is then conducted in the way similar to the classical case. 
We define separation variables {Xi,Pi) by (^. They are canonical if the condition (|^ is 
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satisfied. In the quantum case however, Eqs. (^ are operator equations, hence we have to 
specify the order of operators appearing in A{Xi) and B{Xi). We always assume that the 
position operators precede the momenta. We also assume that all substitutions are done 
from the left. From this it follows that the new momenta Pi are not true observables since 
Pi 7^ P- in general. Now we show how the hermitian separation variables can be defined 
and the separated Schrodinger equations obtained for the Calogero-Moser-Gaudin model 
discussed in this paper. 

We consider the N = 2 case and we take the L-operator (|15D as a Lax matrix. We 
observe that the Gaudin algebra for this model is obtained directly from ([16|) by replacing 
the Poisson brackets { , } with the commutators —i[, ]■ From definition (|^) we find a 
quantum generating function for commuting integrals of motion 



t(A) 


= P'- 


2g2 


M CQ 0/3 _|_ 


.)(A- 




The full set of conserved quantities reads 






H 


= P'- 


2g2 ' 








Ga 


= (Ss) 


1 

i 1 \qa qai 








Ha 


M < 

= >;- 

/3=1 

I3^a 


2^3"^^ + {S^Sl + 


S^_S'^ 


"^ -J- 9 




Ca — e/3 




"T ^. 



+ E 



M 2p^3" + ^([^-,^+]+-[5-5_ 



a=l 



A — er 



(28) 



2q 



[Sl,S4 



and t(A) is given by (p!?! ). Similarly to the classical case, the first of the integrals of 
motion (^) plays the role of the Hamiltonian of the quantum system and Y^a=i Ha = 0. 
For each a = 1, . . . , M, the integral Ga is a quadratic Casimir operator of the algebra 
S0a{2, 1) generated by S°. Hence each of Ga is a number in any irreducible representation 
of the system. The particular representation is obtained by realising spin variables S", 
a = 1, . . . , M in terms of the canonical variables (x^, Pa), [^a, P/s] = i^a/s, a, P = 1, . . . , M 
as follows: 

'-'3 ~ 'ji^aPa +PaXa), '-'+ = '^P^) '-'_ = ^^^a- (29) 

In the representation (^) the first integrals (pSD take the form (compare Eqs.(pOD) 

d2 m 

H - V^^^.^_Pl. 

Ga = 4' « = 1,...,M, 

lb 



l^Ml^ + 112 I 1 ^ 1 ^, 2 2 2 2, 

4 ^^1 e^ - e/3 ^ 4g ^^^ 

Notice that R is an operator now. We can proceed to separation of variables in the 
Schrodinger equations, defining new canonical coordinates by (|^). As we observed, the 
momenta Pj, i = 1, . . . , M are not hermitian, therefore true separation momenta are still 
to be defined. We use Eqs.(0) to find that 

where 

A= "t/1^^~'°\ , ^ = i,...,m. 

The expression for A{X) is derived by the analysis of the behaviour of A{X) at A = 
Xi, . . . ,Xm and A = 00. From ( PB| ) we learn that each Di is positive. Using the fact 
that [Xi, Pj] = i6ij and also that 

^ 1 nf=i(x.-6.) _ nf=i(A-6j 



we find that 



Thus 



-tA-x,n^-,(x,-x,) nfii(A-x, 



^^A) = ^i?(A) [v^^Y^J^Pl^^ 



D,P, = P}Di, (30) 

because A^(A) = v4(A). Therefore we can define hermitian operators 

n, = yAi^-^. (31) 



The operators (|3lD are canonically conjugate to Xi and play the role of true separation 
momenta. Directly from the definition of the Pi we can derive equations of motion (H), 
which in terms of the IT,- read 



1 r.2 fW.,r^ ^ . (r..^( E^ , 3/16 






liiJD,^ (Xi, . . . ,Xm)- E+Y. ^r^^ + 7^ ^ ^ (^1. • • • ,Xm) = 0, (32) 



2 = 1,..., M. Here i?, E^ are eigenvalues of the operators H, Ha, a = 1, . . . , M. To see 
that Eqs.(0) are really separation equations for the model we set 



M 



<iI{X^,...,XM) = ^J\qV\l[<iI,{Xi), (33) 
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where V denotes the Vandermonde determinant V = Yli^j{Xi — Xj). Inserting the wave 
function (|33| ) into equations ( |32D and representing IIj by —i-;^, i = 1, • • • , M, we obtain 



1 d' f /77,. \ , .... , ^ / E^ , 3/16 



/ a=l \^i - ^" 



^/CidX^ \ J ^^^ \Xi — ea [Xi — e 

where 

M 

n i^^ - ^c) 



^, = 0, 2 = 1,...,M, (34) 



a 



a=l 



Equations ( P^ ) are then the separated Schrodinger equations for the generahsed Calogero- 
Moser modeL Finally we notice that the separated equations in the general representations 
of the algebras sOq,(2, 1) can be obtained from (0) by replacing 3/16 with the eigenvalues 
ga of Ga which characterise these representations, i.e. 



Q^,) +E^, + Y, (- 



\/CidX^ V / a=l\-^i~^oi {Xi — ta 

6 Conclusion 

In this paper we have described an integrable generalisation of the Calogero-Moser model, 
which is achieved by coupling the M-particle Gaudin system to the Calogero-Moser model. 
The integrability of the model has been shown by using the recently introduced dynamical 
r-matrix |^. We have also shown that the functional Bethe-Ansatz of [§] can be employed 
to separate the variables in models governed by this r-matrix despite its dependence on 
dynamical variables. The situation is similar to the one discussed in ||ll|. We have used 
this separation of variables procedure in the case in which there are two Calogero-Moser 



particles. We have also used the quantum counterpart of this procedure [|T3] to separate 
variables in the Schrodinger equations. 

The model discussed in this paper is simply an example of a system which is governed 
by a dynamical r-matrix and can be solved by separation of variables in the framework 
of the functional Bethe Ansatz. By an analogy to the non-dynamical r-matrix case, one 
can expect that there are several other models which are governed by dynamical r-matrices 
of the type discussed here, and hence can be solved in a described way. This observation 
opens up new possibilities for constructing and solving integrable models as well as for 
getting a deeper insight into the nature of the models governed by dynamical r-matrices. 
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